1-D Simulation of the Electron Density Distribution in a Novel Nonvolatile Resistive 

Random Access Memory Device 
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The operation of a novel nonvolatile memory device based on a conductive ferroelectric/non- 
ferroelectric thin film multilayer stack is simulated numerically. The simulation involves the self- 
consistent steady state solution of Poisson's equation and the transport equation for electrons as- 
suming a Drift-Diffusion transport mechanism. Special emphasis is put on the screening of the 
spontaneous polarization by conduction electrons as a function of the applied voltage. Depending 
on the orientation of the polarization in the ferroelectric layer, a high and a low resistive state are 
found giving rise to a hysteretic I-V characteristic. The -Rhigh to -Ri ow ratio ranging from > 50% to 
several orders of magnitude is calculated as a function of the dopant content. 

Keywords: Ferroelectric, screening, conductivity, simulation, nonvolatile memory, FRAM, resistive RAM 
(RRAM). 
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Nonvolatile memories based on ferroelectric capacitor 
structures (FRAM) are considered as a potential candi- 
date to replace the current Flash memory generationiS^. 
However, the ongoing trend to miniaturization may level 
the technological efforts in bringing FRAM into produc- 
tion, since charge based devices require a minimum area 
for a save information readout. Various concepts of re- 
sistive memories (RRAM) e.g. based on magnetic tun- 
nel junctions (MRAM)A, phase change materials (PCM, 
Ovonics)£*&, ferroelectric tunnel junctions*, ferroelec- 
tric diodes 8 - 9 - 10 or charge trapping/ detrapping effectsAi, 
which do not show this disadvantage, are currently un- 
der development. The key issue of all approaches is to 
overcome the limitations of today's Flash memories con- 
cerning writing speed, cycling endurance and large pro- 
gramming voltages. 



This contribution reports on numerical simulation 
studies of a memory device consisting of a planar con- 
ductive ferroelectric/ non-ferroelectric multilayer stack, 
which combines the ferroelectric information storage 
principle with a low resistive information readout. In 
the following, it will be referred to as Ferro Resistive 
Random Access Memory, FRRAM. The storage mecha- 
nism is based on the establishment of potential well or 
a potential barrier for electrons, which can be controlled 
by polarization reversal in the ferroelectric layer. This 
causes an enrichment or depletion of the majority charge 
carrier electrons at the ferroelectric/non- ferroelectric in- 
terface, whereby the structure is electrically character- 
ized by a high and a low resistive state. In contrast to 
a ferroelectric diode concept, the proposed device con- 
sists of two highly conductive thin film layers sandwiched 
between two low ohmic contacts. 



Fig. 1(a) shows the capacitor-like memory device struc- 
ture consisting of a ferroelectric and a non-ferroelectric 
layer in contact with metal electrodes. To elucidate the 
operation principle of the device, we first regard the fer- 
roelectric and the non-ferroelectric layer to be ideal insu- 
lators. In a refined model approach, a high concentration 
of conduction electrons will be assumed to be present in 
both layers. 
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FIG. 1: Simplified operation principle of the proposed mem- 
ory device consisting of (a) an ideal ferroelectric/ insulator 
stack, (b) space charge distribution p in the multilayer, (c) 
the resulting inner electric field E and (d) inner potential ip 
for different orientations of the spontaneous polarization. De- 
pending on the orientation of Ps in the ferroelectric layer, a 
potential barrier or a potential well for electrons is formed. 
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A. Ideal Ferroelectric/ Insulator Multilayer 



In the absence of mobile charge carrier in the multi- 
layer, the space charge distribution p{x) for short circuit 
conditions is given by 



p(x) 




x = 
x = l\ 
x = h + h 
else 



(1) 



whereby pf roe denotes the influenced charge on the metal 
electrode and ps is the spontaneous polarization charge of 
the ferroelectric as shown in Fig. 1(b). The relationship 
between space charge and the inner electric field E or the 
electric potential <p are given by Maxwell's first equation 
and the material equation 



dE _ p 
dx e 

and by Poisson's equation 

d 2 (p p 
dx 2 e ' 



(2) 



(3) 



where e is the permittivity of the respective material. 
From the space charge distribution, the inner electric 
field and the inner potential are estimated (see Fig. 1 (c)- 
(d)). It is found that the shape of the inner potential 
is strongly affected by the orientation of the polariza- 
tion. For Ps TT &x we predict the formation of a potential 
barrier, whereas a potential well is formed for Pg fj, e x . 
Here, e x denotes the normal vector in rr-direction. Both 
states can be addressed by application of an external volt- 
age of respective polarity so that the voltage drop across 
the ferroeletric layer exceeds the coercive voltage of the 
ferroelectric material. The present state of the poten- 
tial barrier can easily be sensed by measuring the con- 
ductivity when conduction electrons e.g. by doping are 
introduced into the structure. The latter is expected to 
cause a screening of the polarization charge. Electrons 
will then be attracted by the positive polarization charge 
and repulsed from the negative polarization charge. The 
question, whether or not the shape of the potential esti- 
mated for an insulation system (see Fig. 1(d)) will signif- 
icantly change, if electrons and donor atoms are present, 
will be tackled in the next subsection. We will focus on 
the redistribution of electrons in an n-type semiconduc- 
tor multilayer structure as a function of the orientation of 
the polarization parallel and anti-parallel to an electronic 
current, which is driven by an external voltage. 



B. Conducting Ferroelectric/ Semiconductor 
Multilayer 

If conduction electrons are present in the structure due 
to shallow donor impurities, Eq. (1) has to be replaced 
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FIG. 2: Electron density distribution as a function of the 
applied voltage: I, case of enrichment (potential well). 




+1V 



20 25 30 

x (nm) 

FIG. 3: Electron density distribution as a function of the 
applied voltage: II, case of depletion (potential barrier). 



with 



p(x) 




x = 
x — l± 
x = h + h 
else 



(4) 



Here, en is the elementary charge, n denotes the concen- 
tration of conduction electrons and is the concentra- 
tion of ionized donor atoms. Following a Drift-Diffusion 
approach, the redistribution of electrons under the inner 
electrical field and a concentration gradient is given by 12 



dn 

' dx 



H n nE, 



(5) 



where D n is the diffusion coefficient of the electrons, \x n 
the electron mobility and j n the current density. Under 
steady state conditions, the transport equation simplifies 
to 



dn 

17 



djn 

dx 



0. 



(6) 
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Low ohmic contact properties of the electrode inter- 
faces were approximated by neutral contacts. Then, the 
boundary conditions for electrons are given by 



n(x = 0) = n(x = h+ h) = 



(7) 



The potential difference between the system bound- 
aries is determined by the applied external voltage. For 
the sake of simplicity, a homojunction between the ferro- 
electric and the semiconducting layer and identical elec- 
trodes on both sides are assumed 



ip(x =l% + l 2 ) ~ (f{x — 0) = V e , 



(8) 



Under conservation of charge in the system, the redistri- 
bution of electrons can be calculated as a function of the 
external voltage and the orientation of the spontaneous 
polarization. 



may be suitable for the proposed memory concept. We 
regard complex oxides (Perovskites) 3 , II- VI mixed crys- 
tals (ZnCdS)Ifi, polymers (PVDF-TFE) 13 , but also liquid 
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III. Results 

Numerical simulation techniques are applied to deter- 
mine the self-consistent solution of Eqs. (2)-(8). The elec- 
tron profile for — IV < V cx tcrn < IV for charge carrier en- 
richment is shown in Fig. 2. The case of depletion of elec- 
trons as a function of the applied voltage is illustrated in 
Fig. 3. Both ferroelectric and semiconducting layer have 
a donor concentration of JVd = 10 19 cm -3 and the spon- 
taneous polarization is set to |Pg| = 10/iC/cm 2 . Note: 
A switching of the ferroelectric polarization has not yet 
been considered. For Pg > 0, a strong enrichment of 
electrons leads to a highly conductive state. In case of 
Pg < 0, the depletion of electrons at the ferroelectric/ 
semiconductor interface correlates with a poorly conduc- 
tive state. As shown in the simulation, even if a large 
number of electrons is present in the ferroelectric layer, 
the polarization charge is not perfectly screened, but still 
significantly affects the concentration of electrons in the 
system. Although the inner potential reveals a certain 
distortion of the triangular shape as shown in Fig. Id, 
which is caused by a band bending, the formation of a 
potential barrier or a potential well can still be observed. 

An estimation of the I-V curve from the electron dis- 
tribution under consideration of a switching of the fer- 
roelectric polarization at an external voltage exceeding 
0.5V and an electron mobility of /x n = lcm 2 /Vs is illus- 
trated in Fig. 4. A hysteretic resistance curve is observed, 
whereby the high or the low resistive state can be read 
out by a voltage below the threshold voltage. The inset 
illustrates the scalability of the resisitivity under varia- 
tion of the donor content. Further adjustment can be 
achieved e.g. by variation of the layer thicknesses, the 
value of the spontaneous polarization (different material 
systems) or the contact properties of the electrode inter- 
faces. Even materials with a small spontaneous polar- 
ization, which can not be used in FRAM applications, 



FIG. 4: Estimated I-V characteristic for a donor concentra- 
tion of JVq = 10 19 cm -3 . Inset: Scalability of the resisitivity 
and .Rhigh versus R\ ovi ratio as a function of the dopant con- 
centration. 



crystals as possible candidates for the realization of the 
FRRAM device. 



IV. Conclusion 

We have demonstrated by numerical simulation studies 
the resistive storage properties of a novel memory device 
built from a conducting ferroelectric semiconductor bi- 
layer structure. A high and a low resistive state of the 
memory can be addressed by switching the polarization 
in the ferroelectric layer. Special emphasis was put on the 
interaction between conduction electrons and the ferro- 
electric leading to a redistribution of the majority charge 
carrier and a partial screening of spontaneous polariza- 
tion. An estimation of the I-V characteristic as well as 
the influence of the dopant concentration shows a wide 
scalability of the resistance. This illustrates the poten- 
tial of the proposed device as an alternative concept for 
future high integrated non-volatile random access mem- 
ories ranging from Si-based technology to polymer appli- 
cations. 
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